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Abstract. We present results from an investigation of the clustering evolution of field galaxies between a redshift of z ~ 1 
and the present epoch. The current analysis relies on a sample of ~ 3600 galaxies from the Calar Alto Deep Imaging Survey 
(CADIS). Its multicolor classification and redshift determination is reliable up to I — 23 mag . The redshift distribution extends 
to z ~ 1.1, with formal errors of a z ~ 0.02. Thus the amplitude of the three-dimensional correlation function can be estimated 
by means of the projected correlation function w(r p ). The validity of the deprojection was tested on the Las Campanas Redshift 
Survey (LCRS), which also serves as a "local" measurement. We developed a new method to overcome the influence of redshift 
errors on w(r p ). We parametrise the evolution of the clustering strength with redshift by a parameter q, the values of which 
give directly the deviation of the evolution from the global Hubble flow: £(r comov j ng = lh~ Mpc) = £o(l + z) q . From a 
subsample of bright galaxies we find q = -3.44 ± 0.29 (for ft m = 1, f2 A = 0), -2.84 ± 0.30 (for ft m = 0.2, fi A = 0), 
and q = —2.28 ± 0.31 (for fl m = 0.3, J1a = 0.7), that is a significant growth of the clustering strength between 2 = 1 and 
the present epoch. From linear theory of dark matter clustering growth one would only expect q = — 2 for a flat high-density 
model. Moreover, we establish that the measured clustering strength depends on galaxy type: galaxies with early type SEDs 
(Hubble type: E0 to Sbc) are more strongly clustered at redshifts z <; 0.2 than later types. The evolution of the amplitude of 
the two-point correlation function for these "old" galaxies is much slower (q — —0.85 ± 0.82 for fl m = 0.3, J1a = 0.7). 
Since the evolution of the clustering of bright and early type galaxies seems to converge to the same value in the local universe, 
we conclude that the apparent strong evolution of clustering among all bright galaxies is dominated by the effect that weakly 
clustered starburst galaxies which are common at high redshifts z ~ 1.0 have dimmed considerably since then. Thus the true 
clustering of massive galaxies is better followed by the early types. This provides both a natural explanation for the seemingly 
conflicting results of previous studies and accords with the absence of "faint blue galaxies" in the local universe. 
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1. Introduction 

At the present epoch galaxies are highly clustered and the uni- 
verse seems to be homogeneous only on the very largest scales. 
The question of how the structures we see today, and how the 
galaxies embedded in the dark matter density field have formed 
and developed, is still one of the most challenging ones in the 
field of cosmology. 

First systematic analyses of the distribution of galaxies and 
clusters did not occur before galaxy catalogues with large num- 
bers of objects were drawn up - the first analyses of the cluster- 
ing properties of galaxies were based on the Shane-Wirtanen, 
the Zwicky catalogue, and the catalogue of Abell clusters, and 
the results are outlined in a number of fundamental papers 

j ■ I I r r 

by Peebles (and co-workers) (Peebles 1973; Hauser & Peebles 
1973 | [Peebles & HauseJl974llPeebleJl974llPeebles & Groth 
1 975t IPeebleJ 1 975t) . 

Only a few y e ars later the CfA surv ey was completed 
JPavis et alJ Il982t iDavis & Peebles! Il983l) . and the analysis 
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brought a distribution to light, which was amazingly inhomo- 
geneous - filaments, sheets, walls and large voids emerged, and 
it became clear that the local universe is in fact far from homo- 
geneous. 

The local universe has later been studied using the 
largest catalogues available today, the Las Campanas Redshift 
Survey (Shectman et al. 1996), the Sloan Digital Sky Survey 
JYork et alJ l200(j IStoughton et alJ EooH) . and the 2dFGRS 
icolless etall200ll) . 

It was shown by a number of authors that in the local uni- 
verse Galaxies are biased tracers of mass, i.e. older galax- 
ies are much more strongly clustered than young, starform- 
ing galaxies, b right galaxies are more strongly clustered than 
faint galaxies jPa yis & Gelled 1 197 rl lLahav & Sasla will 9921: 
Santiago A Strauss! M 99* llovino et all ll 993t IfTermit et al. 
1996j [Guzzo et allll997t Lovedav et alJll995l [Willmer et al. 
1998; lNorbergetal.ll2002t IZehavietalJl2002al) . The explo- 
ration of the processes which lead to the different clustering 
of galaxies of different Hubble types can help us to understand 
the interaction between the pure structure growth of the dissi- 



2 



S. Phleps and K. Meisenheimer: Evolution of galaxy clustering 



pationless dark matter component and the development of the 
baryonic matter into stars and galaxies. In principle, the evolu- 
tion of the correlation function can also place constraints on the 
cosmological parameters which determine the geometry and 
dynamics of our universe. Detailed observations of the evolu- 
tion of the clustering strength of different galaxy types have to 
be compared to model predictions (large TV-body simulations 
including starformation and feedback, to disentangle dark mat- 
ter clustering growt h and the evolutio n of the b ias. A first at- 
tempt was made by iKanffmann et ~aTl il999alhh . who carried 
out a semianalytic simulation of galaxy formation and cluster- 
ing in a ACDM cosmology, in which they analysed the cluster- 
ing evolution of galaxies with differing luminosity, color, mor- 
phology and starformation rate. 

For the investigation of structure formation and evolution, 
measurements of the clustering strength extending to redshifts 
of z ft 1 are required. Until recently, there have been only 
few ways to study the large scale structure of the unive rse at 
redshifts z 1, see for example iLe Fevre et alJ lll996l) . and 
ICarlberg et alJ(l2000l) . Not only the shallow depths of most sur- 
veys, but also missing redshift information or too small num- 
ber statistics have limited the possibilities of analysing the 
data with regard to structure formation. In general, two dif- 
ferent types of surveys have to be distinguished - large angle 
surveys, which are limited to relatively bright apparent mag- 
nitudes, and pencil-beam surveys with small, but very deep 
fields. Furthermore, one can distinguish between surveys that 
contain only a small number of galaxies, but with very accu- 
rate redshift information (deduced from spectroscopy), and sur- 
veys that provide huge catalogues of galaxies, but without or 
with rather limited redshift informati on. The Calar Alto Deep 
Imaging Survey (CADIS), see Meisenheimer et al 1 (120021k is 
a deep pencil-beam survey, the output of which at present is a 
catalogue of ~ 6000 classified objects to / < 23 mag . Around 
4000 of these galaxies have reasonably well determined red- 
shifts inferred by means of multicolour methods (a z as 0.02). 
This unique data base provides the possibility to investigate the 
evolution of galaxy clustering from a redshift of z ~ 1.1 to the 
present epoch. 

Usually, structure is described in terms of n-point corre- 
lation functions, the simplest of which is the two-point corre- 
lation function. In practice, computing the three-dimensional 
real-space two-point correlation function requires very accu- 
rate distances. Peculiar velocities as well as redshift errors dis- 
tort the redshift-space relation, and by making the distribution 
more Poisson-like, increase the noise. Different methods have 
been developed to overcome these problems. If no redshifts are 
available at all, it is possible to obtain information about the 
three dimensional distribution of galaxies by deprojecting the 
two-dimensional angular correlation function w(0). If peculiar 
velocities are not negligible, or the data suffer from large red- 
shift errors, one can use the deprojection of the projected cor- 
relation function JDavis & Peeblesll983l) to deduce the cluster- 
ing strength of the three-dimensional distribution. This method 
is used to derive the results of the present paper. 

This paper is structured as follows: The Calar Alto Deep 
Imaging Survey and the data used for the analysis are described 
in Section|2] An introduction into the fundamental principles of 



three- and two-dimensional correlation functions and the de- 
projection method used in this paper is given in Section [3] In 
Section0]we investigate the evolution of the galaxy clustering 
from a redshift of z ~ 1.1 to today, in Section|5]the results are 
discussed. 

2. The Calar Alto Deep Imaging Survey 

The Calar Alto Deep Imaging Survey was established in 1996 
as the extragalactic key project of the Max-Planck Institut fur 
Astronomie. It combines a very deep emission line survey car- 
ried out with an imaging Fabry-Perot interferometer with a 
deep multicolour survey using three broad-band optical to NIR 
filters and up to thirteen medium-band filters when fully com- 
pleted. The combination of different observing strategies facil- 
itates not only the detection of emission line objects but also 
the derivation of photometric spectra of all objects in the fields 
without performing time consuming slit spectroscopy. Details 
of the survey and its calibration will be given in Meisenheimer 
et al. (in preparation). 

The seven CADIS fields measure w 1/30 0° each and are 
located at high Galactic latitudes to avoid dust absorption and 
reddening. In all fields the total flux on the IRAS 100 /im maps 
is less than 2MJy/sr which corresponds to Eb-v < 0.07. 
Therefore we do not have to apply any colour corrections. As a 
second selection criterion the fields should not contain any star 
brighter than w 16 mag in the CADIS R band. In fact the bright- 
est star in the four fields under consideration has an R magni- 
tude of 15.42 mag . Furthermore, the fields are chosen such that 
there should be at least one field with an altitude of at least 45° 
above the horizon of Calar Alto being observable at any time 
throughout the year. Among the CADIS fields three equatorial 
fields allow follow-up observations with the VLT. 

All observations were performed on Calar Alto, Spain, in 
the optical wavelength region with the focal reducers CAFOS 
(Calar Alto Faint Object Spectrograph) at the 2.2 m telescope 
and MOSCA (Multi Object Spectrograph for Calar Alto) at the 
3.5 m telescope. For the NIR observations the Omega Prime 
camera was used. 

In each filter, a set of 5 to 15 individual exposures was 
taken. The images of one set were then bias subtracted, flat- 
fielded and cosmic corrected, and then coadded to one deep 
sumframe. This basic data reduction steps were done with the 
MIDAS software package in combination with the data reduc- 
tion and photometry package MPIAPHOT (developed by H.- 
J. Roser and K. Meisenheimer). 

2. 1. Object detection and classification 

Object search is done on the sumframe o f each filter using 
the s ource extractor software SExtractor (IBertin & Arnoutsl 
1996). The filter-specific object lists are then merged into 
a master list containing all objects exceeding a minimum 
S/N ratio in any of the bands. Photometry is done us- 
ing the program Evaluate, w hich has been developed by 
Mei senheimer and Roser (see [Meisenheimer & Rosen Il986l 
and iRoser & Meisenheimer! Il99ll) . Variations in seeing be- 
tween individual exposures are taken into account, in order to 
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get accurate colours. Because the photometry is performed on 
individual frames rather than sumframes, realistic estimates of 
the photometric errors can be derived straightforwardly. 

The measured counts are translated into physical fluxes out- 
side the terrestrial atmosphere by using a set of "tertiary" spec- 
trophotometric standard stars which were established in the 
CADIS fi elds, and which are cali brated with secondary stan- 
dard stars JOkel 1 99Ct IWalshl 1 9951) in photometric nights. 

From the physical fluxes, magnitudes and colour indices 
(an object's brightness ratio in any two filters, usually given 
in units of magnitudes) can be calcul ated. The CADIS m ag- 
nitude system is d escribed in detail in lWolf et alJ (EoOlc) and 
iFriedet al.N200ll) . 

With a typical seeing of l'/5 a morphological star-galaxy 
separation becomes unreliable at R <; 21 where already many 
galaxies appear compact. Quasars have point-like appearance, 
and thus can not be distinguished from stars by morphology. 
Therefore a classification scheme was developed, which is 
based solely on tem plate spectral energy distributions (SEDs) 
JWolf et al.l 1200 ldbl) . The classification algorithm basically 
compares the observed colours of each object with a colour 
library of known objects. This colour library is assembled from 
observed spectra by synthetic photometry performed on an ac- 
curate representation of the instrumental characteristics used 
by CADIS. 

The spectral lib rary for galaxies is d erived from the mean 
averaged spectra of lKinnev et alJ dl996l) . From these, a grid of 
20100 redshifted spectra was formed covering redshifts from 
z = to z = 2 in steps of Sz — 0.01 and 100 different intrinsic 
SEDs, from old populations to starbursts (SED = 1 corre- 
sponds to an E0 galaxy, whereas SED — 100 is a starburst 
galaxy). 

Using the mi nimum variance estimator (for details see 
IWolf et al.l2001cl) . each object is assigned a type (star - QSO - 
galaxy), a redshift (if it is not classified as star), and an SED. 
The formal errors in this process depend on magnitude and type 
of the object. For the faintest galaxies (/ > 22) they are of the 
order of a z = 0.02, and <jsed = 2, respectively. 

Note that we do not apply any morphological star/galaxy 
separation or use other criteria. The classification is purely 
spectrophotometric. 

Details about the performa nce and reliabilit y of the classifi - 
cation are given in Wolf et al. (20 lei) . and lWolf et al . (2001b). 

Rest-frame luminosities can be estimated directly combin- 
ing the redshift information and flux measurements in the 16 
filters; we do not apply evolutionary corrections. 

Four CADIS fields have been fully analysed so far (for co- 
ordinates see Table [0- We identified 4540 galaxies with / < 
23. Out of these, 3626 are located in the redshift range 0.2 < 
z < 1.07, where we have analysed their clustering properties. 
The number of galaxies per field is given in Table together 
with the number of bright galaxies (Mb < —18 + 5 log h). 

The redshift distribution (see Figure [0 has a maximum at 
z « 0.6. The large variation between adjacent bins, which is 
discernably more than expected from Poisson statistics, reflects 
the large scale structure of the galaxies. With only four fields, 
the accumulation of galaxies at certain redshifts does not can- 
cel out completely. We analyse the clustering properties of the 



Table 1. The coordinates of the four fields investigated so far, 
and the number of galaxies per field, / < 23 and 0.2 < z < 
1.07, also givenis the numberof bright galaxies (Mb < —18 + 
5 log h) in this redshift and apparent magnitude range. 
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galaxies between z = 0.2 and z — 1.07, where the statistics 
are large enough and the redshifts are reliable. 
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Fig. 1. Redshift distribution galaxies with I < 23 in the four 
CADIS fields. The "smoothed" redshift distribution can be de- 
scribed by a Gaussian. The dotted lines indicate the redshift 
range, in which the clustering properties of the galaxies are 
analysed. 

There has never before existed a sample of this depth and 
size, that includes redshift information. However, the redshift 
accuracy is not sufficient to calculate the three-dimensional 
two-point correlation function directly. Nevertheless, as we 
will demonstrate later, we can use the redshifts for our analy- 
sis. First of all it enables us to divide the catalogue into distinct 
redshift bins and we can analyse the clustering of the galaxies 
in each of them up to z f» 1.1. Second, we are able to overcome 
the limitations of the angular correlation function. Instead we 
can use the projected correlation function, which makes direct 
use of redshifts, and is less sensitive to redshift errors than the 
three-dimensional calculation. 

We can further divide our sample into different subsamples 
and investigate the clustering properties of galaxies of different 
SED types, or different absolute restframe B magnitudes. 

3. Three- and two-dimensional correlation 
functions 

The simplest statistic commonly used to describe structure is 
the two-point correlation function £(r). It describes the excess 
probability of finding pairs of galaxies separated by a distance 
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r over the probability expected in a uniform, random distribu- 
tion of galaxy locations. Thus, the number of pairs of galaxies 
dN pair (with one galaxy in the volume element dl^ and the 



other in dV^) is given in the form: 

dN pair = N 2 [l + £(r)]dVi dV 2 . 



(1) 



On intermediate scales the correlation function of galaxies is 
featureless. On very small (kpc) scales a change in the slope 
is expected JZehavi et alJl2002bl) . and on scales 20 to 30Mpc 
a cutoff in the angular correlation function has been reporte d 
jMaddox et alJl990tlCoiiins et alJl993lConnollv et alJ2002l) . 
However, the smooth, featureless behaviour at the intermediate 
scales investigated in the present anal ysis can be para metrised 
by a power law, as first proposed by Totsuii & Kihara ( 1969): 



(2) 



The correlation length ro is a measure for the strength of the 
clustering: the probability of finding a galaxy separated by the 
distance r = ro from another galaxy is twice as large as for a 
random distribution of galaxies. Note that the parametrisation 
by ro is equivalent to determining the correlation amplitude at 
a specific linear scale. 

Since CADIS is a pencil beam survey with a field of view of 
11' x 11', it is of little use to estimate the correlation length ro 
from the measurement, because it is well outside the measured 
range of distances. Instead we will deduce the amplitude of 
£(r) at a comoving separation of r = lh^ 1 Mpc. 

Throughout this work correlation func tions are d etermined 
by using the estimator invented by Landv & Szalav ( 1993): 



(DD) - 2(RD) + (RR) 
(RR) 



(3) 



where (DD) is the total number of pairs of galaxies in a given 
radial bin. (RR) is the histogram of the distances between pairs 
of randomly distributed data points that exhibit the same area 
on the sky, and suffer from the same selection effects as the real 
data, thus acting as a catalogue window, which is representing 
the geometrical properties of the survey. (RD) is the histogram 
of the distances between real and random data. Each of the his- 
tograms is normalised to the total number of pairs counted. 

Rand om error s of £ can be calculated from Poisson statis- 
tics (Baush 1993): 



^Poisson — 



[DD] 



(4) 



where [DD] is the unnormalised histogram of the numb er of 
pairs o f real galaxies in logarithmic distance bins. As Bauah 
( 1996) claims, this may underestimate the errors, but they have 
found that the 2a errors computed in this way are compa- 
rable to_fhe_£rrors obtained with bootstrap resampling meth- 
ods tearrow et alJ l*1984). As a test we calculated bootstrap er- 
rors for one cosmological model and one single redshift bin, 
and found the same result. Therefore throughout this paper 
the errors of all correlation functions are calculated following 

— 2(Tp 2 Sson . 



The real space two-point correlation function can only be 
calculated directly, if redshift information is available with very 
high precision, and if peculiar velocities are negligible. This 
is definitely not the case for the multicolour data of CADIS. 
However, it is possible to derive the parameters by inverting 
two-dimensional distributions. 



3. 1. The angular correlation function 

The angular correlation function w(6) is related to the three 
dimensional correlation function £(r) by Limber's equation 
(Limberll954l) . Limber's equation only holds for local samples 
with z <C 1, For deep samples, which cover a large range of 
redshifts, the expansion of the universe, curvature effects, and 
the possible evolution of the clustering have to be included. 
The redshift-dependent version of Limber's equation h as been 
derived bv lGroth & Peeblesldl977h . and Phillipp seTal] ( 1 1 978t) . 
The general distribution of redshifts (dN/ dz) has to be known 
or calculated from the selection function, which assembles all 
selection effects due to observation and data reduction. 

If £(r) is assumed to be a power law (equation the 
correlation length ro of the three-dimensional distribution can 
be calculated. 

The disadvantage of this method is that the inversion is 
highly dependent on the redshift selection function assumed for 
the survey, which is not a direct observable. Any method which 
makes direct use of the measured redshifts, like the projected 
correlation function, gives much more robust results. 

3.2. The projected correlation function 

If the errors in the redshift measurement are not excessively 
large, but still too large to facilitate a direct computation, one 
can use the projected correlat ion function to derive th e correla- 
tion length ro and the slope 7 fc>avis&Peeblesll 983). It makes 
use of the redshifts, but is much less sensitive to redshift errors 
than the three-dimensional correlation function. 

For small angles r 2 = , 
function is defined as 



7r 2 . The projected correlation 



Hr P ) = 2 J(r 2 +. 2 ) 1 / 2 



dir 



£(r)(r 2 -r 2 p )- x ' 2 rdr 



(5) 



where 7r is the separation along the line of sight. The projected 
physical distance r p between galaxies i and j (the distance per- 
pendicular to the line of sight) can - for a given cosmology - 
be measured very accurately: 



r v = [d A (i) + d A (j)] tan(%/2) , 



(6) 



where d A is the angular diameter distance and #y is the appar- 
ent separation between galaxy i and galaxy j. 
If £(r) = (r/ro)~ 7 , then equation l|5} yields 



j(r p ) = c r y p ^ , 



(7) 
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where C is a numerical factor, which depends only on the slope 



7: 



c = 



r ((7-i)/2) 
r(7/2) 



(8) 



Thus computing w(r p ) provides a measurement of the param- 
eters of the three-dimensional correlation function, namely tq 
and 7. From the amplitude of w(r p ) we can deduce the am- 
plitude of the three-dimensional correlation function at r = 
l/i- x Mpc: 

Withf(r) = (r/r )—i 



£(r = l/i -1 Mpc) = P 



w(r p ) = Cr\ 7 rJ 



= Cr^f (r = l/i _1 Mpc) 



(9) 



(10) 



The expressions to dlOt refer to physical distances as would 
be measured by an observer at the epoch under consideration. 
To facilitate a direct comparison of the possible changes with 
respect to the global Hubble flow, we compute the amplitude 
for a comoving distance r com = \hr l Mpc. With 



£(r com = lft^Mpc) - (r (l + z))T > 

which we will call £iMpc in the following, we find 

w{r p )(l + Z y 



ClMpc 



Cr. 



1-7 



(11) 



(12) 



with the numerical factor C from equation (|8j. 

The errors are calculated following Gaussian error propagation: 



^1 



2 ( d£lMpc 



(13) 



In our case the error of 7 is negligible, as explained in the fol- 
lowing section. The first term completely dominates the error 
of £iMpc> thus the approximation cr\ os ^ 1M = a\ og A is valid 
within a fraction of the statistical error. 

The mean galaxy density is determined from the observed 
galaxy counts in each field, whi ch does not ne cessarily repre- 
sent the the true density (Groth & Peebles 1977). The estimator 
will be on average biased low with respect to the true correla- 
tion by a constant amount 



1 



w tIUB (r p )d 2 Sid 2 S 2 , 



(14) 



where u>t rne (r p ) is the true projected correlation function and 
S is the physical area corresponding to the solid angle of the 
field at the redshift under consideration. If the true correlation 
function is given by equation dl0> . the measurement yields 



i(r P ) 



c r y p 



1-7 



Crl[r x -i -l/{CrZ)} 



(15) 
(16) 



The true amplitude CYq is not known, but X/(CrJ) can be 
estimated by performing a Monte Carlo integration (where we 



use the mean of the pair counts (RR) at a projected distance r p 
of the four fields): 



E [{RR) 



.1-71 



c? 



E(RR) 



(17) 



The true value of CYq can be estimated by fitting equation d!6t 
to the data, taking the value of X/(CVq) from equation HH . 
This value, multiplied by the fitted amplitude CYq, yields the 
integral constraint X. This constant will then be added to the 
measured amplitude A, and from the so corrected amplitude 
the value of £im P c will be inferred by means of equation JT2i . 

3.2.1 . Estimation of the projected correlation 
function 



Following Davis & Peebles one can calculate the pro- 

jected correlation function from 



w{r p ) 



£(r p ,7r) dvr 



(18) 



Since the three-dimensional two-point correlation function has 
the form of a power law, it converges rapidly to zero with in- 
creasing pair separation. Therefore the integration limits do not 
have to be ±00, they only have to be large enough to include 
all correlated pairs. Since the direct observable is the redshift 
z instead of the physical separation along the line of sigh t, w e 
make use of a coordinate transformation dLe Fevre et all 19 96): 



w{r p ) 



where 



+5z 



dn 



Sz az 



(19) 



(20) 



£ = ^ [(! + *) M + 

+(1 - a m - n A ){i + zf + n A 



is the increment in physical dista nce. We calculate th e distance 
along the line of sight following K avser et alJ Jl997l) . 

In order to determine w(r p ) we use the two-dimensional 
correla tion amplitude estimated according to lLandv & Szalavl 
dl993l) . Cesti(?>) in the following: 



Cesti(^p) 



(DD(r p )) 5z - 2(RD(r p )) Sz + (RR(r p )) Sz 
(RR(r p )) Sz 



,(21) 



where the subscript Sz implies that projected distances between 
pairs of galaxies are only counted if their separation in redshift 
space (Az) is less than \5z\. To derive w(r p ), Cesti( r p) has to 
be multiplied with the "effective depth" Arn in which galaxies 
are taken into account: 



w(r p ) = Ccsti(?>) ■ Ary 

r f +5Z dn a 

= Cesti ' / -T- dz 



(22) 



In practice, when dealing with large dz, one has to cope with 
a survey selection function of some kind or another between 
—5 z and +Sz, and not with a top-hat function (of probability 
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unity to find a galaxy within the borders of the survey and zero 
otherwise). That is the survey selection function has to be taken 
into account when calculating Arii . The varying probability of 
finding a second galaxy separated by a redshift Az from a ran- 
domly chosen galaxy can be included in the calculation by mul- 
tiplying the integrand in equation i22\ with the smoothed red- 
shift distribution N z 1 dN/dz, normalised to unity at its max- 
imum (N z is the normalisation constant). Then the integral in 
equation Mil is modified to 



(23) 



r +Sz dTr 


' 1 


dN~ 


Lsz dz 




dz 



Ar|| = / — 



With this correction for the selection function, the "depth" con- 
verges to a fixed value and does not grow infinitely, even if the 
integration limits cover more than the total depth of the survey. 
When comparing our results with those of other authors who 
used the projected correl ation function to deri v e the c lustering 
streng th (see for example lLe Fevre et alJl996l o nCarlberg et alJ 
2000) it has to be kept in mind that they did not take the the in- 
fluence of the survey selection function into account. 

The choice of the integration limit 5z in equation dl 81 or 
equation d23i is somewhat arbitrary. To capture the bulk of the 
correlation signal, the integration limits should exceed the red- 
shift separation corresponding to the correlation length tq, i.e. 
be larger than the redshift corresponding to the pairwise ve- 
locity dispersion, and of course they have to be larger than the 
uncertainty in the redshift determination. 

These requirements set lower limits to the value of Sz. But 
how large a 5z should be chosen? Very large values of 5z might 
more completely integrate the correlation signal, but they do so 
at the cost of considerably increased noise, for two reasons: 
First, the larger the separation of two galaxies along the line 
of sight, the more meaningless (in terms of true distance) the 
projected separation perpendicular to the line of sight becomes. 
Second, in the extreme case that a pair of galaxies is separated 
by physical distances s ^> tq along the line of sight, it is most 
likely not correlated at all, since the correlation function de- 
creases very fast with distance. Nevertheless, such pairs can 
display a very small projected separation and would therefore 
be regarded as strongly correlated. 

In order to find the appropriate integration limits for 
our sample, we used the Las Campanas Redshift Survey, in 
the following L CRS. The LCRS is described in detail by 
IShectman etaD |l996). The survey has a median redshift of 
(z) <~ 0.1, and therefore can be regarded as "local"; the mean 
error in redshift is a z i=s 2.24 • 10~ 4 , that is a cz — 67km s _1 . 

We calculated the projected correlation function of the six 
LCRS stripes for increasing integration limits Sz, fitted the am- 



plitudes at 



0.5ft 1 Mpc (the fit was done in the range 



0.07 < r p < 2ft,- 1 Mpc (-1.15 < logr p < 0.3) to make 
sure we fit where the signal to noise is high) and then calcu- 
lated the weighted mean of the six stripes. Also the choice of 



As can be seen in Figure |2] the estimated amplitude of the 
projected correlation function rises very steeply with increas- 
ing integration limits, and approaches a plateau when pecu- 
liar velocities become unimportant and the undistorted correla- 
tion signal is sampled. The maximum value is reached around 
Sz k, 0.01, corresponding to 3000 kms -1 . Obviously the am- 
plitude of the projected correlation function reaches its asymp- 
totic value at the point where the integration limits have about 
the same size as the typical velocity dispersion in clusters 
(Av ~ 2.36 • a v « 2500km s" 1 ). 

Adding errors in the redshift measurement basically lead to 
increasing noise in the correlation signal. If the clustering in 
redshift space is more and more washed out (the redshift dis- 
tribution becomes more and more Poisson-like), the amplitude 
decreases, especially at small scales. 

To prove this assumption and to estimate the size of the 
effect, we assigned artificial errors to the measured redshifts of 
the galaxies in the LCRS catalogue. The errors were randomly 
drawn out of a Gaussian error distribution: 



z — z + Az , 

where Az is randomly drawn from the distribution 



p(Az) 



1 



271-er 



exp 



1 (Az 

2 1^7 



(24) 



(25) 



We calculated the estimate of the projected correlation function 
for the modified samples, one time with a z = 0.007, and one 
time with a z = 0.017 1 for increasing integration limits Sz. 
For both cases the fitted amplitude of w(r p = 500ft _1 kpc) for 
increasing integration limits is also shown in Figure|3 



- Real, measured redshift 
Artificial errors: <j z =0.007 

- Artificial errors: o 7 =0.017 



-f- 



Fig. 2. The influence of redshift measurement errors on the pro- 
jected two-point correlation function. w(r p ) for increasing in- 
tegration limits is calculated for <r z = 0.007 (dotted line), and 
for a z = 0.017 (dashed line). The errors of the fits are not plot- 
ted here to avoid confusion. The marks indicate the integration 



r p = 500ft 1 kpc allows the most direct comparison with the limits which have been chosen for the calculation of w(r p ) 
CADIS data, where we fit the amplitude at r p w 316ft _1 kpc 
(at the mean redshift of the survey (z) this corresponds to a co- ' representing the core of the distribution of the CADIS redshift er- 
moving separation of « 505h^ 1 kpc). Thus, differences in 7 rors as estimated by the comparison of multicolor redshifts with spec- 
have little influence on the comparison. troscopic redshifts for a sample of 160 galaxies 
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As expected, redshift errors require larger integration lim- 
its for a stable estimation. For the artificial redshift errors with 
a z = 0.007, it can be located at around Sz ~ 0.025, for 
the calculation with a z = 0.017, the maximum is reached at 
Sz sa 0.05, that is in both cases at Sz « 3<r 2 . From these 
simulations we conclude that one starts to sample the corre- 
lation signal when the integration limits are larger than the 
Full Width at Half Maximum of the redshift error distribution 
(Az = 2v2 In 2 • a z ~ 2.3a z ). Thus we conclude that cal- 
culate the projected correlation function for the CADIS data, 
the appropriate choice of the integration limits appears to be 
Sz = ±0.05. 

Note that the estimate for the amplitude is systematically 
lowered in the presence of significant redshift errors. For er- 
rors of a z = 0.017, we find that the maximum amplitude is 
a factor 1 .2 lower than in the case of the unchanged data. For 
a differential comparison between various redshift intervals a 
diminished amplitude is in any case no problem, provided one 
uses the LCRS in its modified form (with artificial redshift er- 
rors simulated exaclty as those found in the CADIS survey) as 
local reference. 



4. The evolution of galaxy clustering from z = l 

4. 1. Clustering properties of all galaxies 

We calculated the projected correlation function for the CADIS 
data in three different redshift bins of similar size, 0.2 < z < 
0.5, 0.5 < z < 0.75, 0.75 < z < 1.07. 

A random catalogue consisting of 30000 randomly dis- 
tributed "galaxies" was generated for each CADIS field with 
the same properties as the real data. The surroundings of bright 
stars in our fields were masked out. The same mask is applied 
to the random catalogue. 

The calculation was carried out for a flat high-density 
model (f2 m = 1, Sl^ = 0), a hyperbolic low-density model 
(f2 m = 0.2, f^A = 0), and a flat low-density model with non- 
zero comological constant (O m = 0.3, f^A = 0.7). 

For the CADIS data the projected correlation function 
w(r p ) was fitted over the range 0.02 < r p < 0.5 h~ 1 Mpc 
(-1.7 < logr p < -0.3), whereas the LCRS data was fitted in 
the range 0.07 ^ r p < 2hr x Mpc (-1.15 < logr p < 0.3), 
see FigureFJ] The fitted amplitudes A at r p = 316ft. -1 kpc for 
CADIS and at r p = 500ft _1 kpc for the LCRS, respectively, 
and the amplitude of the three-dimensional correlation func- 
tion at r com = 1 Mpc derived from them, are listed in Table |2] 
(we have left 7 as a free parameter). 

Although we correct for the missing variance on large scales 
by adding the estimated integral constraint X to the measured 
amplitudes, cosmic variance may still play a role. However, as 
can be seen from Figure 0] the errors of the weighted mean of 
the correlation function (if at all) only slightly underestimate 
the varying clustering strength measured in each single field. 
In Figure0]we show the amplitudes A of the projected correla- 
tion function at r p = 0.316ft _1 Mpc measured in each individ- 
ual field (no integral constraints are added yet), for the bright 
galaxies in the lowest redshift bin, for Q, m = 0.3, Ci\ = 0.7. 




log r p /h 1 Mpc 



f! m = 0.3, n A = o.7 




log r p /h 1 Mpc 



Fig. 3. Measured projected correlation function in three redshift 
bins, for the flat standard model (top), an open model (middle), 
and for a flat universe with a non-zero cosmological constant 
(below). Note that the binning in the highest redshift bin is dif- 
ferent. Note also that the calculation was carried out in proper 
coordinates as measured by an observer at the particular red- 
shift under consideration. The lines are the best fits of the data 
in the indicated range. 



The solid line is the amplitude measured from the weighted 
mean, the dotted lines are the la levels. The spread of the mea- 
sured values around the mean is reasonably small, and within 
their errors they lie within the error of the weighted mean. 

Although the projected distances at which the amplitudes 
are fitted are still close to the largest separations measured, the 



8 



S. Phleps and K. Meisenheimer: Evolution of galaxy clustering 



Q„ = 0.3, n.=o.7 



n m =i. n A -o 
n m -o.2, n A =o 



Fig. 4. The amplitudes of the projected correlation function at 
r p = 0.316/i _1 Mpc estimated in each individual field, for the 
bright galaxies in the lowest redshift bin, for f2 m = 0.3, 51a = 
0.7. The solid line is the amplitude measured from the weighted 
mean, the dotted lines are the ler. levels 



measured amplitude A and the slope 7 are independent, the 
error ct 7 is negligible, and therefore the second term in equation 
Jl 3i is negligible. Then 



dlog glMpc 

d\ogw(r p ) 



"'log w(r p ) 



From equation dl 2i we find 31 ° s ^ 1Mpc 



8 log w{r p ) 



1, SO 



(26) 



(27) 



The values of the integral constraint, which are then added 
to the measured amplitudes of the projected correlation func- 
tion, are listed in Table [2] The amplitudes of the three- 
dimensional correlation function at r com = 1 Mpc are derived 
from the corrected amplitudes of w(r p ). 

4.2. Clustering evolution of bright galaxies 

It is well known that in the local universe bright galaxi es are 
much more strongly clustered than faint ones jNorberg et alJ 
120021) . as expected if bright (massive) galaxies form in the 
high density peaks of the underlying dark matter density field 
(Kaiser 1984). For a reliable determination of clustering evolu- 
tion galaxies with similar absolute luminosities have to be com- 
pared. Based on the photometry between 400 and 1250 nm the 
absolute £?-band luminosities can be determined directly with- 
out any X-correction uncertainty. At redshifts up to z ~ 0.3 
CADIS is dominated by faint (Mb < -18 + 5 log h) galax- 
ies (see Figure [5ji, which supress the correlation signal and 
distort the measurement of the clustering evolution. A con- 
sistent comparison of galaxy clustering at different redshifts 
can be achieved, if we calculate the correlation function in 
the lowest two redshift bins only for galaxies brighter than 
Mb = — 18 + 5 log h. For this purpose, the lowest redshift bin 
has to be confined to z > 0.25 instead of z = 0.2, because 
there are only very few bright galaxies at redshifts smaller 
than 0.3. At z = 1.07 the magnitude limit is slightly fainter 
(M B £-18.6 + 5 log h). 



Fig. 5. Redshift distribution of faint (Mb > -18 + 5 log h) 
and bright (Mb < — 18 + 5\ogh) galaxies, for the different 
cosmologies under consideration in this work. 

The amplitudes of w(r p ) and £iMpc are given in Table[3] 
We parameterise the evolution of the clustering strength 
with redshift by a parameter q, which gives directly the devia- 
tion of the evolution from the global Hubble flow: 



(28) 



thus the parameter q can be deduced from a straight line fit in 
a log£-log(l + z) plot (see Figure |6j. q = —2 is expected 
in case of linear clustering growth. For the flat high-density 
model we find q = —3.44 ± 0.29, for the open model we find 
—2.84 ± 0.30, and for the model with non-zero cosmological 
constant q = -2.28 ±0.31. 

The data point of the LCRS is included in the fit for q. Its 
error is extremely small compared to the CADIS data, and thus 
it dominates the fit. If we exclude the LCRS from the fit, we 
get q = 3.86 ± 0.80, q = 3.52 ± 0.84, and q = 3.12 ± 0.90 
(cosmologies in the same order as above). The slopes are sys- 
tematically steeper, although still the same within their errors. 
However, the reduced \ 2 of the fits including the LCRS is 
smaller than that without. We regard this as a corroboration 
that the connection to the local universe indeed improves the 
determination of q. 

The dependency of the parameter q on the cosmology cho- 
sen for the calculation is as follows: The estimate of w(r p ) as 



calculated from equation ( 1221 is the product of Cesti('p) ( as 
defined in equation (12 1> ). and the integral given in equation 
(1231 . The physical distances r p are stretched (or compressed, 
respectively) corresponding to the cosmology, that means the 
individual histograms (DD), (RR) and (RD) are shifted along 
the r p axis. However, the estimator (, e sti{r p ) is unaffected, be- 
cause all histograms are shifted in the same way. The integral 
Ar|| of course changes with cosmology, thus the dependence 
of the amplitude of the projected correlation function on cos- 
mology is due to the integration along the line of sight. Since 
we fit always at the same physical scale (r p = 316/i _1 kpc 
for CADIS), the values of £iMpc we deduce from w(r p ) show 
the same dependency on cosmology. The differences in Arii 
are increasing with increasing redshift. At our highest redshift 



bin ({z) = 0.91) we expect 



Arm (fi m =0.3,fi A =0.7) 



Ar,|(fi„ 



=1,£2a 



1.6, hence 
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Table 2. The measured amplitude of the projected correlation function in different redshift intervals, at r p w 316/i _1 kpc for 
CADIS and at r p w 500/i _1 kpc for the LCRS (first lines), respectively, the integral constraint X and £iMpc as calculated from 
the I-corrected amplitudes of w(r p ), for different world models. Numbers in italic indicate that this values are inferred from the 
LCRS. 



Model 


z interval 


A 


7 


X 


ClMpc 


Om = 1-0 

n A = 


0.04< z <0.16 


7i.92±}™ 


1.92±0.03 


1.89 


14.06±nH 


0.2 < z < 0.5 




1.83 ± 0.09 


5.53 


4.23±!i1| 


0.5 < z < 0.75 


9.16±t65 


1.89 ± 0.14 


3.20 


3.26±;;:2^ 


0.75 < z < 1.07 


3.88±o;gg 


2.03 ±0.18 


1.06 


l-82±^ 


ft™ = 0.2 
ft A = o 


0.04< z<0.1 


73.12±if 7 


1.93±0.03 


1.88 


/•/./«;• 


0.2 < z < 0.5 


2i.44±i;|? 


1.80 ± 0.09 


6.66 


5.39±B:S 


0.5 < z < 0.75 


12.20±^;?? 


1.92 ± 0.14 


3.69 


4.20±^ 


0.75 < z < 1.07 


5.01±i;l^ 


1.97 ±0.22 


1.30 


2.3i±;j;s 


Om = 0.3 

ft A = 0.7 


0.04< z<0.16 




1.91±0.03 


1.89 


15.49±IM 


0.2<z< 0.5 


23.40±S:S 


1.85 ±0.10 


6.30 


5.52±S;I| 


0.5 < z < 0.75 


15.72±|:^ 


1.97 ±0.13 


4.30 


5.31± >2 


0.75 < z < 1.07 


6.i9±f:6l 


2.11 ± 0.25 


1.33 


2-81± :72 



Table 3. The measured amplitude of the projected correlation function at r p w 316/i _1 kpc for the bright (Mb < — 18 ± 5 log ft.) 
galaxies, the integral constraint X, and £iMpc, calculated from the X-corrected A, for different world models. Note that for the 
bright sample the calculation for the first redshift bin was started at z = 0.25 instead of z — 0.2, because there are only very few 
bright galaxies at redshifts smaller than 0.3. 



Model 


z interval 


A 


7 


X 


£lMpc 


ft A - 


0.25 < z < 0.5 


27.38±!;!£ 


1.98 ±0.14 


8.26 


6.80±i;^ 


0.5 < z < 0.75 


11.14±S:SS 


2.02 ±0.19 


4.12 


4.05±&8* 


0.75 < z < 1.07 


3.88± ;b8 


2.03 ±0.18 


1.06 


i.82±8:SS 


ft m - 0.2 
ft A - 


0.25 < z < 0.5 


31.50±i;^ 


1.85 ±0.17 


10.24 


8.02±tss 


0.5 < z < 0.75 


15.42±g:S 


2.02 ±0.16 


4.66 


5.34±^ 


0.75 < z < 1.07 


5.oi±;:m 


1.97 ±0.22 


1.30 


2.3i±g:3 


ft™ - 0.3 
0\ - 0.7 


0.25 < z < 0.5 


39.46±S# 


1.83 ±0.18 


12.74 


io.oi±; : ^ 


0.5 < 2 < 0.75 


22.21±|;^ 


2.00 ±0.16 


5.85 


7.23±i;^ 


0.75 < 2 < 1.07 


6.19±!:^ 


2.11 ±0.25 


1.33 


2.8i±i:^ 



we expect (and find) the amplitudes of the correlation func- 
tion to differ by the same factor. The integral constraint only 
depends on the measured correlation function (see equations 
J15I to H7\ . It is, though estimated with large errors in the 
highest redshift bin, of the order of 25% of the measured am- 
plitude at 316/i _1 kpc in all cosmologies, thus this additional 
offset does not change the ratio of the amplitudes in differ- 
ent cosmologies. If we consider the "local" measurement not 
to change with cosmology, than we expect the parameter q to 
change (due to the change in £iMpc at the highest redshift bin) 
by An — -, — . '°f 1 ; 6 . . = 0.93 between the matter-dominated 

J ^ log 1.91— log 1.1 

flat cosmology and the flat cosmology with non-zero cosmo- 
logical constant. We find Aq — 1.16. The slightly larger differ- 
ence measured by us is caused by the bin at (z) — 0.625 (see 
Figure|6] 

4.3. Comparison with other work 

Before analysing the SED type dependent evolution of galaxy 
clustering, we compare our results with previous attempts to 



study the clustering evolution with redshift. In the literature 
there are essentially only two investigations of the evolution of 
galaxy clusterin g, that can be compare d with the present work: 
one analysis bv iLe Fevre et alJ (11996) which has been carried 
out in the framework of the Canad a France R edshift Survey (in 
the following CFRS), and one bv lCarlberg et alJ J2000h . done 
on the CNOC sample (Canadian Network for Observational 
Cosmology). 

ILe Fevre et alJ ll 19961) used the projected correlation func- 
tion to investigate the spatial clustering of 591 galaxies be- 
tween 0.2 < z ^ 1.1, in five CFR S fie lds (for a descriptio n 
of the survey see Lilly et alJ i 19951) and ISchade et ail dl995h . 
respectively) of approximately the same size as our CADIS 
fields. The objects are primarily located in three parallel strips 
for each of the five fields, within which almost 100% spec- 
troscopic sampling was obtained, separated by regions where 
few spectroscopic observations were carried out. The galaxies 
with spectroscopic redshifts have / < 22. 5 mag . They computed 
the projected correlation function in three redshift bins between 
0.2 < z < 0.5, 0.5 < z < 0.75, and 0.75 < z < 1.0, with in- 
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Q m =i, n.=o 




q = -3.44+/-0.29 



0.1 0.2 
log (1+z) 



6 

S 0.5 



n m =o.2, n A =o 



q = -2.84+/-0.30 



0.2 
log (1+z) 



6 

° 0.5 



n m =o.3, n A =o.7 




q = -2.28+/-0.31 



0.1 0.3 
log (1+z) 



Fig. 6. The evolution of the clustering strength (the amplitude 
of the three-dimensional correlation function at at a comoving 
separation of lh^ 1 Mpc) with redshift. The integral constraint 
X has been taken into account. Filled symbols are only bright 
galaxies, open symbols include the complete sample. The line 
is the fit to the bright galaxies, the first data point (from the 
LCRS) is included in the fit. 

tegration limits of S z = ±0.0075. F or the connection to z = 
they took values of rn (z — 0) from lLovedav et al.l dl995) and 
lHudon&LiiivlJl996l) . 

Figure0shows the amplitude of the three dimensional cor- 
relation function at r com = 1ft, -1 Mpc, deduced from the pro- 
jected correlation function, in comparison with our own data. 
Since they did not apply corrections for the missing variance, 
we calculated £im P c without adding I to the measured ampli- 



tudes of w(r p ). For direct comparison, we have to multiply our 
measured amplitudes of the projected correlation function of 
the complete sample by jpf^jfgg] = 1-2 in order to allow 
for the large errors of the CADIS multicolour redshifts (see 
Section|3j- With this correction, the CFRS data points are con- 



ro=5.0 
7=1.7 



n m =i, n.=o 



log (1 + z) 

Fig. 7. The amplitudes of the three dimensional correlation 
function at r com = 1ft -1 Mpc, deduced from the projected 
correlation function of the CFRS data (crosses), in compari- 
son with our own data (open symbols). The amplitudes are cal- 
culated without taking the integral constraint i nto account, to 
allow for a direct comparison with the data of iLe Fevre et alJ 
(1996). The data points with integral constraint correction can 
be seen in Figure[6]for comparison. The filled symbols are our 
data points corrected for the influence of the redshift errors on 
the amplitude of the projected correlation function (shifted by a 
factor of 1.2). The dotted line is the fit of the CFRS data points 
including the value of tq(z = 0) from lLovedav et aDdl995h . It 
seems that they underestimated the errors, since from a sample 
of roughly six times fewer galaxies they derived errors of the 
same size as ours. 

sistent with our own measurement for the sample as a whole. 
However, it seems that they underestimated the errors, since 
from a sample of roughly six times less galaxies they derived 
err ors of the same size as ours. 

ILe Fevre et alJJl996ll claim that if rn(z = 0) = 5/i _1 Mpc, 
< e < 2. The fit of their data points, including a direct con- 
nection to z = 0, that does not take the different properties of 
the samples into account, formally yields q — —3.04 ± 0.21. 
If the connection to z = is disregarded, we find q = 
-1.18 ± 0.63 from their data. 

This exercise shows that the value of q derived from the 
CFRS solely depends on the connection to the present epoch. 
To deduce a reliable result, it is indispensable to take into ac- 
count that the catalogues to be compared have to consist of the 
same mix of Hubble types and luminosities, as we have seen, 
bright galaxies are more strongly clustered. It is not possible to 
estimate the evolution of the clustering strength by comparing 
a sample of intrinsically faint galaxies at redshifts z <; 0.3 with 
the bright galaxies which dominate the local measurements. 
The comparison between different samples has to be carried 
out in a selfconsistent way. The integral constraint, which they 
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did not take into account, is different in different redshift inter- 
vals (see Table|2ji, and thus alters the result. Also the measured 
amplitude of the projected correlation function depends on the 
redshift accuracy and the appropriate choice of the integration 
limits. 

Thus we conclude that their measurement is consistent with 
ours, but they did not treat the local measurement selfconsis- 
tently^ 

ICarlberg et all §000) carried out an analysis of the clus- 
tering evolution on a sample of 2300 bright galaxies from 
the CNOC surv ey. The survey itself is described in detail in 
lYee et all dl99€l) . The galaxies have fc-corrected and evolution- 



k.e 



< 



-20 11 



(H = 
= 0.65. 



compensated R luminosities M\ 
100km s _1 ). The redshift distribution extends to z 
For comparison with z = 0, they selected a comparable sample 
from the LCRS. They also do not apply corrections for the inte- 
gral constraint, but since the survey area of the CNOC is large 
(extending over 1.55 D° in total) X is expected to be negligible 
at the scales the y investigate (0.16 < r„ < S.O/t -1 Mpc). The 
parameter e (see lCirofh & Peeblesil977tlEfstathiou et alJl99lh 
that they use for the parametrisation of the clustering evolution 
can of course be related to the parameter q 



g=-(3 + e-7) 



(29) 



With equation J29l > we can calculate the parameters q which 
correspond to the e values they give: They found q = —2.00 ± 
0.22 for Vt m = 1, O a = 0, q = -1.37 ± 0.18 for Q ro = 0.2, 
Q A = 0, and q = -0.39 ± 0.19 for Q m = 0.2, fi A = 0.8. 

Their values indicate a much slower evolution of the clus- 
tering than our results. The reason for the descrepancy might 
be the different sample selection - the CNOC galaxies have 
been selected to be bright in the R band, whereas the galaxies 
in our sample have bright blue luminosities. As we will show in 
the next section, the clustering growth of blue and red galaxies 
evolves differently, and in fact the q values we derived for our 
subsample of early type gala xies are within the err ors identical 
with the values estimated by Carlber g et alJ J2000l) . 



4.4. The evolution of the correlation function for 
different Hubble types 

Various local surveys have established that the 



red 



blue ones JDavis & Geller 


1 19761 lLahav & Saslawl 1 19921 


Santiago & Strauss 1992 


Lu 


ovino et al. 1993; Hermit et al. 


1996; Guzzo et al. 1997; Willmer et al. 1998). 


lLeFevreetalJJl996) 


separated their sample into red and 



blue galaxies (bluer and redder than the Colem an et alJ 41980) 
Sbc spectral energy distribution), and found that they have 
comparable clustering amplitues at z > 0.5. Here we show 
that this is not the case, at redshifts z ^ 0.2 red galaxies are 
more strongly clustered than blue ones. 

The galaxy library used for the multicolour classification 
resembles regular grids in redshift and SED, thus the Hubble 
type can also be estimated from the observations. This enables 
us to investigate the evolution of the clustering of different pop- 
ulations of galaxies. 



We devided the sample into two SED bins, with SED < 
60, including galaxy types E to Sbc (1355 galaxies), and 
SED > 60, including Sds to Starbursts (23 1 1 galaxies), re- 
spectively. Galaxies with earlier SEDs have smaller redshift 
errors, therefore we primarily concentrate on the investigation 
of the evolution of the large scale structure of the galaxies in 
the SED < 60 sample. 

Figure|8]shows w(r p ) for the different redshift bins, in com- 
parison with the projected correlation function of all galaxies, 
for a flat fl\ = 0.7 model. The early type galaxies show signif- 
icantly stronger clustering than the late types. Table|4]lists the 
amplitudes at r p = 316/i _1 kpc, fitted between ~ 20h~ 1 kpc 
and ~ 500/i -1 kpc, and the integral constraint I. 



f! m = 0.3, C! A =0.7 




— 0.2<z<0.5 
e 0.5<z<0.75 
• 0.75<z<1.07 

1 1 1 I 1 

-2 



log r p /h 1 Mpc 



n m =o.3, n A =o.7 



— > — 0.2<z<0.5 
o 0.5<z<0.75 
0.75<z<1.07 




~ 2 log r p -/h-' Mpc 

Fig. 8. The projected correlation function in three different 
redshift bins; upper panel: early type galaxies in comparison 
with all galaxies (only the fits are shown as grey lines), lower 
panel: late type galaxies in comparison with all galaxies for 

n m = 0.3, n A = 0.7. 



The corresponding amplitudes of the three-dimensional 
correlation function £(r) at a comoving distance of r com = 
Ih^ 1 Mpc, corrected for the missing variance I, are given in 
Table|5] 

Figure [9] shows the amplitudes of the three-dimensional 
correlation function £(r) at a comoving distance of r = 
lh- 1 Mpc, for red galaxies (SED < 60), for Q, m = 0.3, 
Qa — 0.7. For comparison the data points for the whole sample 
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Table 4. The measured amplitude of the projected correlation function at r p = 316h 1 kpc and the integral constraint X, in 
different redshift intervals for two SED bins, for different world models. 



Model 


z interval 


4(S'£D < 60) 


7 


1 


A(SED > 60) 


7 


I 


o in 

Qa ~ 


0.2 < z < 0.5 


46.67 ± g;ii 


1.95 ±0.09 


14.33 


7.80 ± mi 


1.96 ±0.24 


0.14 


0.5 < z < 0.75 


21.33 ± $:S 


2.13 ±0.14 


6.12 


4.97 ± Hi 


2.01 ±0.33 


0.66 


0.75 < z < 1.07 


5.53 ± 3 ; 42 


2.84 ±0.93 


1.01 


4.78 ± i:S§ 


1.98 ±0.22 


0.57 


a« - 0.2 

Q A - 


0.2 < z < 0.5 


59.17 ± fcg 


1.97 ±0.09 


15.59 


9.09 ± 1% 


1.92 ±0.27 


2.13 


0.5 < z < 0.75 


27.97 ± 


2.01 ±0.16 


9.18 


5.89 ± ITo 


2.20 ±0.44 


1.31 


0.75 <z< 1.07 


11.03 ± 


2.69 ±0.80 


1.42 


5.67 ± 


0.52 ±0.86 




Om - 0.3 

Ov - 0.7 


0.2 < z < 0.5 


67.60 ± ?:ta 


2.20 ±0.10 


15.14 


6.62 ± t%l 


2.37 ±0.30 


1.19 


0.5 < z < 0.75 


39.22 ± fcgg 


2.00 ±0.11 


11.12 


10.40 ± ii§ 


2.16 ±0.44 


2.30 


0.75 < z < 1.07 


19.53 ± ftg 


2.16 ±0.64 


3.35 


9.25 ± 


1.89 ±0.28 


1.98 



Table 5. The amplitudes of the three-dimensional correlation function £(r) at a comoving distance of r = 1 Mpc, for red 
(SED < 60) and blue (SED > 60) galaxies; the integral constraint has been taken into account in the calculation. 



Model 


z interval 


CiMpc (S££> < 60) 


6m P c (SBD > 60) 


n m - 1-0 
n A - 


0.2 < z < 0.5 


n.26±i:S§ 


l-46±g:S 


0.5 < z < 0.75 


7.27±s:ii 


i.49±g:is 


0.75 < z < 1.07 


2.34±J:3 




n m - 0.2 
n A - 


0.2 < z < 0.5 


13.77±J;?9 


2.08±g:iJ 


0.5 < z < 0.75 


9.87±i;S§ 


i-89±6:is 


0.75 < z < 1.07 


4.56±?;g^ 




fim - 0.3 
- 0.7 


0.2 < z < 0.5 


i5.n±i:i? 


i.29±ss 


0.5 < z < 0.75 


is.37±i:Si 


3.35±1;?| 


0.75 < z < 1.07 


8.57±§:2i 


4.oi± ( 1 ,:^ 



are also plotted. The amplitudes are corrected for the missing 
variance X. The lines are the fits for the q-parameter (equation 
d). 

The late type galaxies have large redshift errors, and they 
are intrinsically only very weakly clustered. Although their 
number density is high, the large errors and the low correlation 
signal result in a highly uncertain measurement of their cluster- 
ing amplitudes (the measurement is shown for the Q m = 0.3, 
Oa = 0.7 case), which would make a reliable estimate of the 
clustering evolution impossible. Also the late type sample is at 
high redshifts dominated by relatively bright galaxies, whereas 
at lower redshifts the rather low-mass, fainter galaxies, which 
are only very weakly clustered, are currently undergoing peri- 
ods of high starformation rates. A galaxy classified as extreme 
late type at z = 1 may have evolved into an Sa at z = 0, thus we 
are not dealing with a homogeneous sample. A measurement 
of the clustering evolution of late type galaxies is therefore not 
meaningful at this point. 

The LCRS measurement is not included in the fits for the 
early type galaxies. Table [6]lists the parameter q for the early 
type sample, and compares them to the values found for the 
bright sample. The evolution of the amplitude of the correlation 
function is clearly differential, the clustering strength of early 
type galaxies grows much slower. However, it looks as if the 
clustering evolution of the CADIS early type sample converges 
to the same value. In Table |6] we list also the parameter q for 
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1 1 1 1 


III 




= 0.3, Q A = 0.7 








q = -0.34+-028 
q = -0.85+-O.S2 








q = -2.28+-0.31 
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Fig. 9. The evolution of the clustering strength (at a comoving 
separation of lh^ 1 Mpc) with redshift, for early type (SED < 
60). The data for the bright galaxies and the complete sample 
are plotted for comparison. A correction for the integral con- 
straint is included. The first data point (the weighted mean of 
the LCRS sectors) is not included in the fit for the early type 
galaxies. We also show the fit which does take it into account 
(dotted line), and the measurement for the late type sample for 
the Q m = 0.3, ft\ = 0.7 case. Note that the errors of the late 
type galaxies are too large to measure the clustering evolution. 
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Table 6. The parameter q for early type (SED < 60) galaxies, for a fit done without taking the the LCRS point into account, and 
for a fit including it. The evolution of the bright sample is shown for comparison. 

Cosmology I q(SED < 60), without LCRS q{SED < 60), with LCRS gbright 

f2 m = 1.0, fi A = 0.0 -2.60 ±0.91 -1.52 ±0.30 -3.44 ± 0.29 

Q, m = 0.2, 0\ = 0.0 -2.03 ±0.82 -0.74 ± 0.28 -2.84 ± 0.30 

tt m = 0.3, fi A = 0.7 -0.85 ±0.82 -0.34 ± 0.28 -2.28 ± 0.31 



the early type sample which was determined when taking the 
LCRS point into account. 

When the LCRS point (at z — 0.1) is included in the fit, the 
CADIS points still lie well on the fitted line, which indicates 
that the clustering strength of bright and early type galaxies 
indeed seems to converge to the same point at z ^ 0. 

These values of q, -1.47 ± 0.37, -0.71 ± 0.31, -0.34 ± 
0.32 (see Table [6] for the corresponding cosmological mod- 
els) are within their errors identical to the values deduced by 
ICarlberg et all J2000I) from the red CNOC galaxies (they found 
q = -2.00±0.22,<7 = -1.37±0.18, q = -0.39±0.19forfhe 
cosmologies in the same order as above, see previous section). 

5. Discussion 

We now want to discuss the most important results of our 
present analysis in turn: 

(1.) Bright galaxies are more strongly clustered than faint 
galaxies and show significant growth in clustering since z ~ 1. 

At redshifts around z ps 0.3 the CADIS data are dominated 
by rather faint galaxies, which show smaller clustering ampli- 
tudes than brighter ones. On the other hand, the LCRS sample 
consists mainly of intrinsically bright galaxies. In order to be 
able to compare our high redshift data self-consistently with 
the local measurement from the LCRS, we measured the am- 
plitude of the correlation function from a subsample of bright 
galaxies (Mb < —18 + 51og/i). We find a significant growth 
of the clustering strength for this bright subsample: 

For the Q m , fl\ = 0.3,0.7 cosmology we find q = —2.2 
that is ! 1Mpc > z ~"? = 5.5. This might be compared with the the- 

£lMpc(z=l) | 2 . j . 

oretical predictions by Kauffmann et al. ( 1999a). They use N- 
body simulations of dark matter (DM) to derive the clustering 
of DM halos at < z < 5 and semi-analytic models in order 
to assign a galaxy of certain type and luminosity to each DM 
halo. Although many properties of their mock galaxies pop- 
ulation are only in loose accordance with observations, their 
result seems robust: the clustering strength of bright galaxies 
(Mb < —19 + 5 log h) on 2 Mpc scales follows closely that of 
the DM halos they are imbedded out to z ~ 1 (bias parameter 
~ 1). They find a 3.5-fold increase of ^2Mpc between z = 1 
and z = 0. If one takes into account that the increase should 
be higher on 1 Mpc scales and that our limit for "bright galax- 
ies" is 1 magnitude fainter (i.e. it includes less massive halos) 
it seems that the the clustering evolution predicted by ACDM 
models is in agreement with our measurements. 

(2.) At high redshifts, early type galaxies are more strongly 
clustered than the bright galaxies and the clustering strength of 



the early type galaxies evolves much more slowly than that of 
the bright galaxies. 

A plausible explanation for the different evolution of 
the clustering properties of the early typ e galaxies arises in 
the context of biased galaxy formation jBardeen et alj fl986: 
iBrainerd & Villumsenlll994li . The first galaxies are born in a 
highly clustered state, because they form in the bumps and wig- 
gles which are superimposed on the very large-scale density 
enhancements. The next generations of galaxies form later in 
the wings of the large-scale enhancements, and are therefore 
less and less clustered. While the universe expands, the galax- 
ies evolve, age, and eventually merge to form larger, brighter 
galaxies and ellipticals, and generally add to the population of 
earlier type galaxies, while new galaxies form at later times in 
less and less clustered environments. Merging creates galaxies, 
which suddenly" add to the old population. A merger event also 
reduces not only the number of galaxies, but also the number 
of small pair separations in a sample, which reduces the prob- 
ability of finding pairs of galaxies at small distance s - and thus 
supres ses the amplitude of the correlation function. lFried et all 
( 2001) found the density evolution of the early type and the late 
type galaxy population in the CADIS galaxy sample suggestive 
of merging. 

Although the clustering strength of the underlying dark 
matter density field increases with redshift, the biasing de- 
creases. The net effect is a very slowly rising clustering am- 
plitude of the early type (old) galaxies. 

Thus, the measured rate of the clustering growth depends 
on the mixture of galaxy types one observes at different red- 
shifts. 

(3.) The evolution of the clustering amplitude of early and 
bright galaxies converges to the same value at z « 0. 

The convergence of the evolution of the clustering of bright 
galaxies and galaxies with SEDs earlier than Sbc towards the 
same local measurement can also be understood if one takes 
into account that galaxies evolve and hence the population mix 
one observes changes with redshifts. The comoving number 
density of weakly clustered starburst galaxies increases with 
increasing redshifts, whereas the space density of the highly 
clustered very early type (E-Sa) gal axies decreases b y a fac- 
tor of ~ 1.6 from z = to z = 1 JFried et alJl200 lh . In our 
highest redshift bins the clustering signal of the bright sub- 
sample is presumably dominated by a weakly clustered pop- 
ulation of galaxies, which at that time were bright, and blue. 
At z = these weakly clustered blue galaxies have vanished, 
and now the majority of the early type galaxies (which had 
later, bluer SEDs at higher redshifts), are also bright, thus at 
the present epoch there is a large overlap between the bright 
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and the early type samples. Thus we conclude that the apparent 
strong growth of clustering of the bright galaxies is dominated 
by the fading of the unclustered, blue population. 

To understand exactly how galaxy evolution influences the 
measurement of the growth of structure, and the evolution of 
the large scale structure influences the evolution of galaxies, 
we need detailed investigations of the evolution of both the cor- 
relation function and the luminosity function of galaxies with 
different SEDs. 

Larger, wide angle deep surveys have only recently be- 
come available, and one of them is the COMBO 17 survey 
(Classifying Obj ects by Medium-Band Observations in 17 fil- 
ters, fWolf et alJ J2001al) l in some respect the successor of 
CADIS. The complete catalogue will include ~ 40000 galax- 
ies with / < 23, in 1D°, with SED and morphological infor- 
mation. This amazing data base can be used for various inves- 
tigations, using the projected correlation function. The higher 
statistic allows for a more detailed analysis of the evolution of 
the clustering of different galaxie types and their relation to the 
underlying dark matter density field. 

The clustering properties of starburst galaxies at higher 
redshifts can be investigated using the emission line galaxies 
observed by the CADIS emission-line survey using an imag- 
ing Fabry-Perot interferometer (Hippelein et al. 2002). These 
galaxies, which have been detected and classified by their 
emission lines, have redshifts with an accuracy of 120 km 
s _1 - good enough to calculate the three-dimensional corre- 
lation function directly. The special observing technique sam- 
ples galaxies in distinct narrow redshift bins, which allows for 
the investigation of the evolution of the clustering properties 
of emission line galaxies between a redshift of z w 1.4 and 
z w 0.24. 

The observations have to be compared to theoretical pre- 
dictions. We plan to carry out a large cosmological simulation 
including starformation and feedback. The individual galaxies 
will be assigned an SED according to their stellar masses and 
starformation histories, so we can perform synthetic photome- 
try applying our CADIS or COMBO 17 filter set to the syn- 
thetic spectra. Thus we will be able to "observe" the mock 
galaxies at different redshifts, and directly compare the sim- 
ulations to our observations. 
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